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SLIDE  1:  (1.5  minutes) 

Both  the  design  and  the  analysis  of  underwater  acoustic  systems  are 
dependent  on  good  models  of  the  local  acoustic  ambient.  It  Is  well  known  that 
the  sea  surface  generates  acoustic  radiation  that  often  dominates  the  system 
output  noise  field,  and  despite  considerable  effort,  Improvements  In  the  mod¬ 
els  remain  necessary  when  hydrophone  arrays  having  vertical  apertures  are  used 
In  a  wide  range  of  environments.  The  objective  of  the  task  that  1  am  discuss¬ 
ing  Is  the  development  of  a  model  of  the  vertical  directional  spectrum  due  to 
that  searsurface-generated  radiation.  Directional  spectra  has  been  chosen 
because  It  Is  a  readily  visualized  statistic  that  Is  physically  Intuitive  and 
has  proven  useful  In  the  design  and  analysis  of  underwater  acoustic  systems. 

The  present  state  of  the  theory  of  sea  surface  sound  Is  such  that  the 
model  development  Is  necessarily  experimentally  based. ^  While  several  mea¬ 
surements  of  the  directivity  function  have  been  made.  It  has  proven  difficult 
to  make  this  measurement  In  such  a  way  that  the  near  horizontally  directed 
energy  Is  not  significantly  affected  by  long-distance  sources  that  complicate 
the  modeling  process.^  A  theme  In  the  data  processing,  which  I  will  describe. 
Is  that  having  the  total  directional  spectra  determined  by  the  local  sea  sur¬ 
face  condition  allows  one  to  partition  the  measurement  into  more  than  a  single 
source  type  because  the  source  types  distinguish  themselves  by  their  direc¬ 
tional  patterns. 

The  approach  this  project  took  was  to  build  a  database  of  directional 
spectra  measured  In  an  acoustically  Isolated  area  where  the  acoustic  ambient 
Is  dominated  by  the  local  sea  conditions. 
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SLIDE  2:  (1  minute) 

The  Tongue  of  the  Ocean  (TOTO) ,  The  Bahamas,  is  the  location  of  the 
measurement.  It  is  a  relatively  deep  water  basin  that  is  totally  isolated 
from  global  shipping.  This  figure,  with  the  depth  contours  in  fathoms,  shows 
that  the  shallowest  area  is  700  fathoms,  or  1300  meters,  deep.  This  is  deep 
enough  in  wavelengths  not  to  act  as  a  shallow  water  location  but  shallow 
enough  to  have  a  lossy  “bottom- limited”  propagation  characteristic  within  the 
basin  itself.  The  Naval  Underwater  Systems  Center  has  a  presence  in  TOTO 
which  makes  available  environmental  monitoring  locations  and  ready  access  to 
small  vessels  for  servicing  the  acoustic  measurement  system. 

Two  measurement  locations  have  been  used.  The  initial  two -month  deploy¬ 
ment  of  the  measurement  system  was  at  a  near  reef-site  close  to  Fresh  Creek. 
The  location  was  chosen  for  logistic  reasons  to  service  the  initial  wetting 
down  of  the  system.  The  present  location  of  the  system  is  on  the  east  side  of 
the  cul-de-sac.  This  location  is  further  shielded  from  transiting  vessels  in 
TOTO  by  the  land  mass  north  of  the  location.  Ten  months  of  the  planned  one- 
year  deploymenthave  been  completed.  Note  that  the  diameter  of  the  cul-de-sac 
is  a  little  more  than  30  nml. 
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SLIDE  3:  (2  minutes) 

The  measurement  system  objective  was  a  broadband  characterization  of 
the  diffuse  acoustic  ambient  field.  The  Naval  Underwater  Systems  Center 
IR/IED  sponsor  required  that  the  system  be  quite  low  cost.  The  system  that 
you  see  in  this  figure  is  the  resulting  compromise.  The  dilemma  of  the  acous¬ 
tic  antenna  design  is  that  over  much  of  the  angular  range  the  diffuse  ambient 
changes  slowly  with  vertical  angle  and  a  low- resolution  system  is  adequate, 
but  near  the  horizontal  direction,  where  media  refraction  can  abruptly  change 
the  solid-angle  density  function,  a  considerably  higher  resolution  system  is 
necessary.  The  fielded  system  covered  the  two-decade  frequency  band  shown 
with  a  modest  2-  to  4 -wavelength  aperture.  A  parametric  spectral  estimation 
algorithm  that  models  the  media  refraction  was  then  used  to  obtain  the  higher 
angular  resolution  near  the  horizontal  direction. 

The  seven  octaves  are  recorded  sequentially  by  synchronously  recording 
the  eight  14 -bit  hydrophone  signals  from  each  section  on  three  video  cas¬ 
sette  recorders.  Also  recorded  are  the  depth  and  tilt  of  the  vertical 
hydrophone  array.  After  the  21-minute  recording  period  the  system  sleeps 
for  50  hours  before  it  repeats  the  cycle.  Every  45  to  50  days  it  is  neces¬ 
sary  to  retrieve  the  system  to  replenish  the  3  VCR  tapes  and  96  flashlight 
batteries  that  power  the  electronics. 

The  brightly  colored  syntactic  foam  collar  around  the  instrumented 
pressure  vessel  maintains  the  verticality  of  the  system  within  a  small  frac¬ 
tion  of  the  hydrophone  array  angular  resolution. 
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SLIDE  4:  (1  minute) 

The  basic  statistic  measured  is  the  hyd^'ophone  cross -spectral  density 
(CSD)  matrix.  However,  because  the  CSD  matrix  is  not  a  readily  vlsxialized 
quantity,  it  is  useful  to  calculate,  from  the  CSD  matrix,  a  beamformed  output 
that  is  displayed  and  examined  to  deteirmine  the  quality  of  each  measurement. 
This  figure  Illustrates  the  “direction  of  arrival  spectrum”  in  (dB///iPa^/Hz) 
as  a  function  of  frequency  and  elevation  angle.  The  contour  plot  on  the  left 
gives  quantitative  details,  and  the  surface  plot  on  the  right  gives  a  qualita¬ 
tive  view  of  the  entire  surface.  Note  that  zero  degrees  is  looking  straight 
up  at  the  water  surface.  The  figure  is  a  composite  of  all  seven  octaves;  the 
high  frequency  end  of  each  aperture  is  visually  “flagged”  by  the  aliasing  lobe 
that  occurs  at  the  frequency  for  which  the  hydrophones  are  spaced  one -half 
wavelength  apart.  Note  also  that  in  the  horizontal  direction,  i.e.,  at  an 
elevation  angle  of  90  degrees,  a  clear  refractive  notch  is  evident  at  high 
frequencies  despite  the  low  angular  resolution  of  the  measurement  system. 
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SLIDE  5;  (2.5  minutes) 

The  overall  characteristics  of  the  vertical  directional  spectra  as  a  func¬ 
tion  of  the  sea  surface  condition  are  illustrated  here.  These  figures  repre¬ 
sent  first-order  estimates  of  the  directional  spectrum  formed  by  simply  nor¬ 
malizing  the  beamformed  output  by  the  solid-angle  coverage  of  the  major  lobe 
and  then  averaging  the  spectral  levels  over  one -third- octave  bands.  The  ordi¬ 
nate  of  these  figures  is  the  directional  spectrum  in  (dB//pPa^/Hz/steradian) . 
The  two  Independent  variables  are  again  frequency  and  elevation  angle,  with 
zero  degrees  looking  up  at  the  water  surface.  Except  for  the  refractive 
“notch”  area  surrounding  the  horizontal  direction,  these  first-order  estimates 
are  quantitatively  quite  correct  and  were  obtained  in  the  initial  near -reef 
deployment.  Higher  order  estimates  of  the  directional  spectra  will  be  dis¬ 
cussed  shortly  but  are  not  necessary  for  the  present  discussion. 

The  measurement  used  to  describe  the  sea  surface  condition  is  a  nondimen- 
slonal  surface  friction  velocity.^  Vind  speed  measurements  made  10  meters 
above  the  water  surface  were  converted  to  surface  friction  velocity  and 
then  nondlmenslonallzed  by  the  minimum  phase  velocity  of  the  water  surface 
capillary- gravity  wave  field.  A  unity  value  of  the  nondimens ional  friction 
velocity  is  interpreted  as  the  point  of  incipient  whitecapping. 

The  following  aspects  of  the  data  should  be  noted.  First,  with  white- 
caps  present  the  measured  vertical  directional  spectrum  in  the  direction  of 
the  sea  surface  is  a  smooth  surface  with  a  rather  flat  maximum  in  the 
direction  of  the  surface  at  500  Hz .  This  surface  extends  over  the  entire 
frequency  range  of  the  measurement.  Second,  without  whitecaps  present  but 
with  nondimens ional  friction  velocity  greater  than  one-half,  the  shape  of  the 
measured  directional  spectrum  surface  above  its  maximum  value  at  500  Hz 
remains  unchanged  although  the  height  of  the  surface  falls  with  the  third 
power  of  the  friction  velocity.  Below  200  to  AOO  Hz  the  surface  falls  rapidly 
with  a  reduction  in  frequency.  Third,  with  the  nondlmens ional  friction  velo¬ 
city  less  than  one -half,  the  measured  directional  spectrum  is  dominated  by  the 
modal  behavior  of  the  TOTO  channel,  not  the  local  sea  surface;  i.e.,  the  sea 
surface  radiation  is  no  longer  being  measured.  Fourth,  below  100  Hz  the  data 
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are  frequently  contaminated  by  “tones”  radiating  over  a  large  horizontal 
extent  of  the  escarpment.  The  last  observation  is  characteristic  of  the  near¬ 
reef  site  data  only. 
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DIRECTIONAL  SPECTRA  WITH  WHITECAPS  PRESENT 
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SLIDE  6;  (3  minutes) 

There  are  three  reasons  for  choosing  a  parametric  spectral  estimation 
approach.  First,  if  an  appropriate  bases  function  set  exists  then  quite 
efficient  parameterization  of  a  complex  field  variable  statistic  is  pos¬ 
sible.  Second,  an  increase  in  spectral  resolution  over  nonparametric  algo¬ 
rithms  is  possible.  Third,  physical  interpretation  of  data  can  be  enhanced 
if  physical  interpretation  of  the  model  is  possible.  For  all  these  reasons, 
a  model  fit  to  the  data  was  attempted.  The  model  development  is  ad  hoc  and 
the  components  of  the  model  occur  historically  in  the  acoustic  ambient  lit¬ 
erature.  This  illustration  is  an  attempt  at  quickly  summarizing  the  ingre¬ 
dients  of  the  model- fitting  procedure.  The  actual  mathematics  are  very 
straightforward. 

All  three  of  the  components  pictured,  or  at  least  some  modified  version 
of  them,  have  been  used  individually  in  other  studies.*’^  What  is  unique  to 
the  present  approach  is  that  more  than  one  source  type  is  assumed  present 
and  the  data  are  then  used  to  assign  a  strength  to  each  one  individually. 

The  three  components  are  (1)  a  plane  of  statistically  independent  vertical 
dipoles  located  at  the  sea  surface,  (2)  a  volume  of  statistically  indepen¬ 
dent  monopoles  exponentially  distributed  in  depth  below  a  smooth  pressure 
release  surface,  and  (3)  an  isotropic  term  intended  as  a  “catchall”  component 
to  Include  such  things  as  scattering  from  media  nonhomogeneities.  The  first 
two  sources  are  propagated  to  the  receiving  system  using  a  multipath  expansion 
of  a  wave  equation  solution.  This  allows  the  modeled  terms  to  be  given  an 
angular  distribution.  This  was  realized  using  the  “generic  sonar  model” 
software  package.  Depending  on  the  frequency,  either  BLUG  or  MGS  bottom  loss 
values  were  used.  The  effects  of  the  hydrophone  spatial  sampling  and  finite 
aperture  are  Included  in  the  model. 

The  parameters  representing  the  source  area  density  occur  linearly  in 
the  model,  and  the  exponential  depth  scale  occurs  nonlinearly.  Because  the 
volume  and  surface  source  terms  become  indistinguishable  from  each  other  as 
the  depth  scale  goes  to  zero,  an  iterative  parameter  estimation  scheme  was 
used  to  ensure  orthogonality  in  the  parameters.  The  result  is  that  a  volume 
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source  term  becomes  significant  only  if  the  depth  scale  is  of  the  order  of 
a  wavelength  or  greater. 

Note  also  that  if  a  contaminating  plane  wave  is  present  in  the  data 
it  also  is  readily  accounted  for  in  the  fit  to  data  but  is  ignored  in  the 
reconstruction  of  the  directional  spectrvim. 
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PARAMETRIC  MODEL  SUMMARY 
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SLIDE  7:  (1  minute) 

An  important  indication  of  the  quality  of  the  model  is  how  well  it  fits 
the  data.  This  figure  shows  an  example  of  the  model -data  comparison  at  the 
beamformed  output.  In  each  panel  the  abscissa  is  the  elevation  angle  steered 
to  and  the  ordinate  is  the  beamformed  output.  In  the  direction  of  the  sur¬ 
face,  l.e.,  between  0  and  90  degrees,  which  is  the  process  we  wish  to  model, 
the  fit  is  nearly  always  within  1  dB  of  the  data.  The  standard  error  for  the 
“fit  to  data”  as  well  as  the  standard  error  of  parameter  estimates  is  typi¬ 
cally  less  than  5  percent.  The  cross-correlation  coefficient  between  parame¬ 
ter  estimates  is  always  less  than  one- tenth. 

Note  that  in  the  upper  left-hand  panel  there  existed  a  “target”,  i.e.,  a 
plane  wave  contaminating  the  data.  A  plane  wave  component  added  to  the 
model  readily  fit  the  data. 
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MODEL-DATA  COMPARISON 


iki| 


““ 

r: 

iVi 

\WM 

»i 

1^1 

■< 

( 

z 

J. 

i- 

— 

— 

— 

— 

— 

— 

— 

V 

_ 

w 

_ 

_ 

_ 

_ 

1^41 


Him 

INSil 

\m\ 


im 


Mt 


— ijij- - 

tm 

rnmm 

- 1^11 - 


■inBaiiHHunB 


amm 


iiiiiviiiiiiiiiP 


TD  8733 


SLIDE  8:  (3/4  minute) 

This  figure  shows  an  example  of  a  measured  vertical  directional  spectrum 
employing  the  parametric  spectral  estimation  described.  The  measurement  is 
from  the  cul-de-sac  deployment  with  no  whitecaps  present.  As  before,  the 
ordinate  is  the  solid-angle  density  function  illustrated  as  a  function  of  log 
frequency  from  40  to  4000  Hz  and  elevation  angle  from  straight  up  to  straight 
down.  The  direction  of  the  sea  surface  is  on  the  right-hand  side.  The 
refractive  “notch”  with  its  frequency -dependent  depth  is  quite  evident. 
Because  the  notch  depth  increases  with  frequency,  it  is  most  likely  due  to 
TOTO  modal  response  rather  than  scattering. 


L8 


In  this  figure  I  have  plotted  the  dipole  coefficient,  which  is  a  source 
area  density  term,  as  a  function  of  frequency  with  nondimens ional  friction 
velocity  as  a  parameter.  The  following  two  characteristics  are  worth 
noting:  First,  for  nondimens ional  friction  velocity  greater  than  one-half, 

the  source  density  is  inversely  dependent  on  the  first  power  of  frequency 
above  1  kHz  and  independent  of  frequency  below  1  kHz  until  a  lower  frequency 
cutoff  is  reached.  The  source  strength  is  proportional  to  the  third  power  of 
the  friction  velocity  without  whitecaps  and  the  second  power  with  whitecaps. 
The  second  observation  is  that  for  nondimens ional  friction  velocity  less  than 
one-half,  either  the  source  strength  is  falling  faster  than  the  third  power  of 
the  friction  velocity  and  is  below  the  background  level,  or  the  dipole  compo¬ 
nent  is  only  erratically  present. 
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SLIDE  10:  (1  minute) 

I  will  now  summarize.  The  objective  of  this  project  is  to  model  the 
vertical  directional  spectra  due  to  sea  surface  sound  alone.  The  approach 
taken  is  to  accumulate  a  database  obtained  in  an  acoustically  isolated  area  in 
which  the  acoustic  ambient  is  dominated  by  the  local  sea  conditions.  To  do 
this,  a  low-cost  vertical  hydrophone  array  operating  from  40  to  4000  Hz  has 
been  operating  in  the  Tongue  of  the  Ocean,  The  Bahamas,  for  over  a  year.  A 
data  processing  system  has  been  designed  and  implemented  that  features  a  par¬ 
ticularly  efficient  parameterization  of  the  measured  directional  spectra.  The 
physically  based  parametric  model  used  in  the  spatial  spectral  estimation  pro¬ 
cess  has  given  some  anecdotal  evidence  of  the  character  of  the  sea  surface 
sound  source  structure.  A  statistically  based  examination  of  the  data  is 
underway.  The  measurement  period  will  be  completed  in  July  of  this  year. 

So  far,  140  trials  have  been  recorded. 
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Low-Cost  Broadband  Vertical  Hydrophone  Array  And  Data  Recording  System 
Designed  And  Installed  In  An  Isolated  Area. 
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